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Abstract 

The energy spectrum of the ^^C nucleus with {J'^,T) = (O+jO) and (2"'",0) is 
investigated in the framework of the multicluster dynamical model by using a deep 
Qa-potential with forbidden states in the S and D waves. A very high sensitivity of 
the compact ground and first excited 2f states energy levels to the description of the 
two-body forbidden states wave functions has been estabilished. It is shown also that 
the chosen method of orthogonalizing pseudopotentials yields convergent results for 
the energies of the excited (0^,0) and (03^,0) states of the ^^C nucleus with a well 
developed cluster like structure. 
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1 Introduction 



A systematic study of the ^^C nucleus structure has been successfully performed on the basis 
of the microscopic three- alpha models [1, 2, 3, 4, 5]. These models reproduce the general 
features of the low- lying ^^C spectrum. In the calculations using the rcsounating group 
method (RGM) and the orthogonality condition model (OCM) the ground and first excited 
2| states energies of the ^^C nucleus are overbound by 7-10 MeV, although the energy 
positions of the 02", 22 , 3j" and 1" resonance states are reproduced fairly well. In the recent 
work [6], which used the RGM method, the low-lying energy spectrum of this nucleus was 
described with a good accuracy. This work has shown also that the astrophysical significant 
0^ state of the carbon is a genuine three-alpha resonance in the continuum. This fact is 
important, since almost all the carbon is produced via cupture reaction a +^ Be — C* + 7, 
which goes through O2 resonance of the ^^C nucleus. 

Although the microscopic models are useful tool for the study of the structure of light 
nuclei, but they can not yield the wave functions in a convenient form for future application 
in nuclear reactions with these nuclei, such as above cupture process. Therefore, it is impor- 
tant to develope macroscopic models. These models should provide main properties of the 
structure of the nucleus to be investigated. 

There are two alternative local potential models, describing the alpha-alpha data. The 
first potential suggested by Ali and Bodmer (AB) [7] involves a strong repulsive core. The 
second deep potential model was proposed in [8] by Buck, Friedrich and Wheatley (BFW), 
which has three forbidden states in the lowest S- and D- waves. Thus, these alternative 
models differ each from other in describing the Pauli repulsion part of the alpha-alpha 
interaction. As a result, the local on-shell equivalent potential models give still different 
wave functions for the ^Be ground state: while the BFW potential yields a nodal behavior, 
the AB potential does not describe this microscopically substantiated property. 

Macroscopic three-alpha models also are beheved to describe well the low-lying states of 
the ^^C nucleus due-to large binding energy of a particle. However, in sharp contrast to the 
microscopic models, the local AB potential gave a very weak binding with E = —0.60 [9] 
for the g.s. energy Eg^p = —7.27 [10] of the ^^C nucleus as a srtong off-shell effect of the 
repulsive core presented in this potential model. On the other hand, when using the BWF 
alpha-alpha potential, one has the problem of elimination of forbidden states from relative 
motion wave function of the system. However, the construction of the many body projector 
is a very complicated procedure even for the case of the three body system. Therefore, 
various approaches to the full many-body projector are used. 

The approach used in the work [11] is based on the strong truncation of the functional 
space of relative motion. In spite of the truncation of the relative motion space, the calculated 
value overbound the ground state binding energy of the ^^C nucleus by about 5 MeV. The 
inclusion of full relative motion space yields a stronger ovcrbinding by about 15 MeV [12] with 
neglecting the Coulomb forces, which is estimated by 5-6 MeV. To overcome the three-alpha 
overbinding problem the authors of the work [12] proposed the Finite Pauli repulsion model 
(FPRM) for the aa interaction. However, this model allows admixtures of the forbidden 
states into the physical relative motion space. 

On the other hand, from above mentioned 3a macroscopic calculations with orthogonality 
conditions, one can not see, how accurate the full 3-body projector is approximated by chosen 
orthogonalization method. Since the antisymmetrization has strong dynamical effects, the 
3a problem has to be investigated in a model, which allows to examine the convergence of 
the orthogonalizing procedure. In other words, the dynamical effects of the Pauli principle 
has to be described properly. 



In the present work we investigate the Set-problem by using the BFW potential in the 
framework of the multicluster dynamical model [13] for light nuclei based on a high accuracy 
variational method which has been used successfully in many structure calculations of various 
atomic and nuclear systems [14, 15, 16, 17, 18]. For the elimination of forbidden states we 
use the method of orthogonalizing pseudopotentials (OPP) which allows to work in the full 
relative motion space. The main feature of our work is the description of the convergence of 
the orthogonalizing procedure for the ground and lowest excited states of the ^^C nucleus. 
When using the OPP method one can examine the convergence of the three-body energy 
as a function of the projecting constant A [16, 19]. We check also the the convergense in 
respect to the description of the two-body forbidden states fixed by chosen aa potential. We 
show that the convergence of the orthogonalizing procedure has different character for the 
compact shell-model like bound states and for the resonance states with a well developed 
cluster like structure. 

In Section 2 we describe the model for the investigation of the 3a-system structure. The 
numerical results are presented in Section 3. Discussion and conclusion are given in Section 
4 and Section 5, respectively. 

2 Model 

For the description of the lowest states of the ^^C nucleus we use the multicluster 3q; model 
with a local aa-interaction potential of Buck, Friedrich and Wheatley [8] 

V{r) — Voexp{—rir^) , (1) 

with Vo=-122.6225 MeV, r] — 0.22 fm^^ [11]. This potential describes well the experimental 
phase shifts of the aa-scattering 6l{E) with L = 0, 2, 4 up to 40 MeV energy. The Coulomb 
interaction potential in our calculations is taken in the form 

Vcouiir) = Ae^erfibr)/r, (2) 

where b=0.75 fm~^, which corresponds to the a particle charge distribution being the Gaus- 
sian form with a width of 1/6. We use a value h'^/rria = 10.4469 MeV fm^ in our calculations. 

The potential (1), with including the Coulomb interaction in the form (2) has three 
nonphysical bound states forbidden by the Pauli principle in each aa-subsystem, with the 
energies E{Ot) = -72.6249149 MeV , E(O^) = -25.6174 MeV, E{2+) = -21.9991037 MeV. 
The of state corresponds to the shell configuration s^, while the 0^ and 2+ states correspond 
to the s^p^. These forbidden states in our three-body calculations are eliminated by using 
the method of orthogonalising pseudo-potentials (OPP). This method was developed in the 
work [20] and successfully employed for the investigations of the structure of nuclei with 
^ = 6 and 3N system [15, 16, 19]. 

The variational method on a Gaussian basis used for the calculations of the spectrum of 
the ground and excited states of the ^^C nucleus has been given in details in Ref.[14]. A high 
accuracy of the method has been demonstrated in a number of works [14, 15, 16, 17, 18]. 
A general formalism of the method and analytical expressions of the matrix elements of the 
overlapping integral, the kinetic energy operator, central, spin-orbital and tensor interaction 
potentials for the system of three identical fermions with a spin value 1/2 have been given 
in Ref.[14]. A corresponding part of these matrix elements can be used for the variational 
calculations of the Sa-system. Therefore, we give here only main formulas for the expression 
of the wave function of the Set-system in series of symmetrized Gaussian basis functions. 



The three-body Schrodinger equation is solved by using a pscudopotential of the form 

Vi{r) = V,{r) + K,„,,,(r) + (3) 

/ 

where A/ is the projecting constant, V^P is the projecting operator to the /-wave forbidden 
state in the two-body subsystem {j + k), {i,j,k) = (1,2,3), and their cychc permutations. 
We note that the method of OPP for the ehmination of forbidden states uses the first term 
of the expansion for the full three-body projector [2] 

p=Y.h-jzhPj+ i: hP.Pk----. (4) 

i=l i,j=l i,j,k=l 

where 

P^ = J:W\ (5) 

/ 

which corresponds exclusively to the two-chister Pauli forces. We neglect three-cluster 
(triple) Pauli forces when using this approach. The method of OPP allows us to obtain 
the solution of the Schrodinger equation at large values of A/ which is orthogonal to the 
two-body forbidden states. 

Thus, we use next three-body pseudo-hamiltonian in our Sa-cluster variational calcula- 
tions: 

H^Ho + Vi + V2 + V3, (6) 

where Hq is the kinetic energy operator of the three a- particles. The wave function of the 
3a- system is expanded in the series of symmetrized Gaussian functions [14] : 



*f^ = E4'V;-, (7) 

13 

where ipf^j = ipyj{l; 2, 3) + ip^j{2; 3, 1) + ^^^(3; 1, 2), 



(pjj{k; I, m) = Njxly[exp{-axjxl - (ii^yD^i^ {xk, m) (8) 

Here {k; I, m) = (1; 2, 3), (2; 3, 1), (3; 1, 2), 7 = (A, J, M) = (70, J, M);xk, Vk are the normal- 
ized Jacobi coordinates in the A;-set: 

Xk = -rm) = r n^m] 

yk = ^( — ^ rk) = T^ Pk, (9) 

Nj is a normalizing multiplier. The nonlinear variational parameters axj^Pij are chosen as 
the nodes of the Chebyshev grid: 

2 J — 1 71" 

axj = aotg{^^-),j = 1,2,...A^a, 

(3ij ^Potgi^^^),]^ 1,2,. ..Ni, (10) 

where ao and /?o are scale parameters for each (A/) partial component of the complete wave 
function. When we use the Chebyshev grid the basis frequensies axj,Pij cover larger and 



larger intervals around the scale parameters as the numbers Nx and Ni increase. This 
allows us to take into account both short-range and long-range correlations of particles. The 
extraordinary flexibility of the many-particle Gaussian basis makes it possible to describe 
three- and four-particle configurations that are formed in the ground state and excited states 
of multicluster systems, and which exhibit an extremely high degree of clustering [14]. 

The angular part of the Gaussian basis (8) is factorized into the angular component and 
the internal wave functions of the a-particles and has the form: 

^xt'i^kjk) = {Yx{ik)<S}Yi{mjM<f>{i)mm (11) 

Here the orbital momenta A and I are conjugate to the Jacobi coordinates Xk and yk, respec- 
tively. 

We note that the use of the symmetrized Gaussian basis allows us to take into account 
contributions actually of all partial waves to the binding energy of the system due to over- 
lapping of various three-body channels with the same total orbital momentum L ( and total 
spin S in the case of the 3N-system). On the other hand, the solution of the Schrodinger 
equation should be symmetric over the permutation of any two a particles. This requirement 
is satisfied automatically when using a symmetrized basis for any number of basis functions, 
whereas in the case of nonsymmetric basis this requirement is fuUfilled approximately and 
depends on the precision of the expansion. 

One of the main advantages of the variational method on a Gaussian basis is the pos- 
sibility of fully analytical calculation of moments of interaction potentials of some special 
forms. For example, to calculate the integral 

POO 

Gia,j3,n)= / x^expi—ax^ — I3x)dx,a> Q,n>Q, (12) 

which arises when Yukawa type potentials are used, we suggest more convenient and effective 
recurent formulas (compare with Ref.[17]): 



Gia,m = ^eM'^)erH^) (13) 

G(a,fl,l) = ^-^G(a,tlM (14) 
la la 

n 3 
G(a, /3, n + 1) = —G{a, /3, n - 1) - ^G{a, P,n),n>0 (15) 
la la 

And for the calculation of moments of the nonpoint Coulomb interaction potential (2) 

/■oo 

Gerf{oi, P,n) = / x'^exp{—ax'^)erf{Px)dx,a>0,n>0, (16) 



we obtained also next recurent formulas: 

G„,(a, P, n) = A » - 2) + ^^J^^, (18) 

These formulas allow us to essentially economize computer time when we tabulate the values 
of G(a, /9, n) and Gerfic^, n) for the fixed values of a and f3. We note also that the zero-th 
moment Ger/(a,/3, 0) can be estimated numerically. However, in our calculations we need 
moments with only odd values of n. 



3 Numerical results 



First we note that a very large symmetrized Gaussian basis is used in our calculations. For 
the calculations of the energy levels with {J^,T) = (0^,0) we used a basis containing 80 
functions in the each of the three-body channels (A, /) = {(0, 0); (2, 2); (4, 4)}. We calculated 
the levels with {J'", T) = (2+, 0) on the basis containing 437 Gaussians chosen in main three- 
body channels (A, /) = {(0, 2); (2, 0); (2, 2); (2, 4); (4, 2); (4, 4)}. The results indicate that a 
further extension of the basis does not have a remarkable influence on the accuracy of the 
expansion. 

The projector on the f-wave forbidden state in each two-body subsystem has the form: 
with the forbidden state function 

Here ro is the "projector radius", and A^^^ is the normalizing multiplier: 

^(/)_2/+7/4 _ 2//2 2 N 

7rV4[(2A + l)!!]i/2' «m-T/^zro^;. 

In order to check the behaviour of the 3«-system energy when improving the accuracy 
of the expansion of the two-body aa forbidden states wave functions in a series of Gaussian 
functions we choose several sets of the Gaussian approximations. Corresponding approximate 
values of the deep forbidden state energies of the aa-system are shown in the Table 1. In 
Set 1 the forbidden state wave function of the O^'-level of the «a-system is approximated 
via N=l Gaussian, and wave functions of the forbidden Oj- and 2+-levels are approximated 
via N=2 Gaussians, et al. The values of the forbidden two-body states energies are given in 
corresponding squares. By comrarison of these numbers with the corresponding exact values 
of the two-body 05*", O2, 2+-forbidden states energies one can conclude about the quality of the 
approximation for the given number N. We note that the wave function of the aa forbidden 
02'-statc contains a node, therefore it can not be expressed by a single Gaussian. 

The spectrum of the energy levels of the ^^C nucleus with (J'^,T) = (0"'",0) and (2"'",0) 
for the several variants of the Gaussian expansion for the forbidden state wave functions 
are shown in Tables 2 and 3. For the sake of convenience we use an identical value of the 
projecting constant A/ for all forbidden states. In all three tables the symbol "< P >" 
denotes the total admixture of the forbidden states to the energy of the 3Q;-system: 

< p >=< I A(A + A + A) I *s > ■ 

In Table 4 we give the experimental values of the lowest energy levels of the nucleus 
and our theoretical estimation comparing with results from the literature. The results of 
the Ref.[ll] were obtained by using the hyperspherical garmonics method and the BFW 
act-potential. Microscopical calculations in the Ref.[6] were performed with an effective MN 
nucleon-nucleon potential. 

From Tables 2 and 3 one can see that we have a reasonable saturation of the energy levels 
of the ^^C nucleus when increasing the projecting constant A to the infinity for a given set 
of approximation of the two-body forbidden states. However, beginning from the value of 



the projecting constant A > 10*^ MeV the ground and first excited 2^ states energies of the 
^^C nucleus begin to display a high sensitivity to the description of the two-body forbidden 
states wave functions. When improving a quality of the approximation of the act-forbidden 
states in these cases wc go to the weakly bound system . 

For the other excited states (02^, Ojj", 22")of the ^^C nucleus we have a quite reasonable 
convergence of the energy with respect to the description of the aa-forbidden states. 

We note also that admixture of the forbidden states in the ground and first excited 2f 
states of the ^^C nucleus decreases slowly as increasing the projecting constant A to the 
infinity. In the case of the 3A^ system at A = 10^ MeV the value of the forbidden states 
admixture in the "^iJ ground state energy was approximately 10~^ MeV when using the 
Moscow NN potential with forbidden states [15, 16]. 



4 Discussion 

First of all we discuss the reason of the strong dependence of the energies of the ground and 
excited 2^" states of the ^^C nucleus on the expansion accuracy of the two-body forbidden 
states wave functions. We know from the experiment that these states are bound strongly 
and have consequently a well developed compact shell-model structure [10]. Thus, in the 
states with (J'^,T) = (01^,0) and (J'^,T) = {2i,0) three a-clusters with a large probability 
overlaps simultaneously each with other. In other words, three-cluster Pauli forces should 
play an important role in the description of these states due to the requirement of the Pauli 
principle. However, the method of OPP used here for the elimination of forbidden states in 
the full three-body system takes into account really two-body Pauli effects. 

It is known that the two-body projectors Pi and Pj do not commute each with other 
and are not mutually orthogonal [2]. We suggest that the overlapping of these two-body 
projectors is very strong for the three-alpha system, and hence the energies of the compact 
ground and excited 2f states of the ^^C nucleus display a high sensitivity to the description 
of the forbidden states beginning from A ~ 10^ MeV where the next terms in the expansion 
(11) begin to play a role. And these neglected terms present a mathematical expression of 
the Pauh triple forces. 

We note that the authors of the Ref.[21] also came to the conclusion that three-body Pauh 
forces should be taken into account when investigating the structure of the 3a-system. They 
made this conclusion on the basis of the calculations of the ^^C nucleus ground state energy 
in the framework of the Alt-Grassberg-Sandhas equation method by using aa-potentials 
obtained by the RGM and modified by the OCM. 

On the other hand, as we noted in the previous section, the energies of the excited levels 
with {.r,T) = (0^,0) and (J^,T) = (0^,0) change slowly from Set 1 (A^ < 2) to Set 4 
(A^ = 7). Moreover, at A 10^ — 10^ MeV one can see a good saturation of the energy of 
these levels with respect to values of the projecting constant A. Thus, for the energies of 
these levels we can write (see Table 2) : Etheori^t) ~ 1-551 MeV and Etheori^t) — 4.055 
MeV. It is important to note that the difference of these values Atheor — 2.50 MeV is in a 
good agreement with it's experimental value A^xp = 2.65 MeV. [10]. And this agreement 
can be seen for all sets 1-4 at large values A ~ 10^ — 10^ MeV. Thus, the excited levels of the 
^^C nucleus with {J'^,T) = (02", 0) and {J'",T) = (03,0) can be described qualitatively in 
the framework of the method of OPP in contradistinction to the ground and first 2f excited 
states. 



5 Conclusions 



The energy spectrum of the ^^C nucleus with (J^,r) = (0+, 0) and (J'',T) = (2+,0) was 
calculated in the framework of the multicluster dynamical model using the ««-potential 
of Buck, Friedrich and Wheatley with forbidden states in the S and D waves. For the 
elimination of forbidden states in the full three-body system the method of orthogonalising 
pseudo-potentials (OPP) has been used. 

It was shown that the energies of the compact (O^", 0) and (2^'", 0) states of the ^^C nucleus 
with a shell-model like structure display a very high sensitivity to the description of the two- 
body forbidden states. Therefore, for the energies of these strongly bound states, the chosen 
model does not give convergent results. We suggest that this situation is due-to neglecting 
of three-body Pauli forces, which should play an important role in the description of the 
structure of these strongly bound states. 

We obtained also that in our model the resonance states (02", 0) and (03',0) with a well 
developed cluster-like structure can be described quite well. In contrast to the ground and 
first excited states, the orthogonalization procedure yields convergent results. However, the 
energy position of the astrophysically significant O2" resonance EexpiS^t) — 0.3796 MeV is 
overestimated by about 1.2 MeV. 
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Tabic 1: Sets of Gaussian approximations for the forbidden states wave functions of the aa 

system and eorrcsponding cncrg}' vahics in McV 



0^ 


0+ 


2+ 


exact E = -72.624915 


exact E = -25.6174 


exact E = -21.999104 


N =1: E = -72.5445 


N = 2: E = -25.106 


N = 2: E = -21.676 


N = 3: E = -72.6126 


N = 3: E = -25.5558 


N = 3: E = -21.8837 


N^A-.E^ -72.6233 


AT = 4 = -25.6111 


N^A-.E^ -21.9576 


N ^7 -.E^ -72.624905 


N^7 :E^ -25.6173 


N ^7 :E^ -21.999098 



Table 2: The energy spectrum of the ^"^C nucleus with {.r,T) = (0+,0) in MeV for the 
several sets of Gaussian approximations of the two-cluster forbidden states at several values 
of the projecting constant 





A (MeV) 


10 


102 


10^ 


10^ 


10^ 


10^ 


10^ 


10« 


< 2 


E, 


-210.80 


-45.194 


-20.665 


-18.999 


-14.882 


-12.85 


-12.56 


-12.44 




E2 


-150.72 


-20.037 


-0.677 


+0.274 


1.210 


1.284 


1.39 


1.42 




E3 


-109.37 


-15.878 


+ 1.239 


1.306 


+2.98 


3.446 


3.99 


4.01 




<P> 


29.83 


30.36 


0.847 


1.133 


1.73 


0.300 


6.E-2 


5.E-2 


N^3 


El 


-210.71 


-43.298 


-19.776 


-16.492 


-6.387 


-3.666 


-3.35 


-3.311 




E2 


-150.05 


-17.183 


-0.411 


-0.021 


1.238 


1.331 


1.42 


1.463 




E3 


-109.53 


-15.652 


+1.344 


1.380 


+2.553 


3.361 


4.07 


4.046 




<P> 


29.93 


28.78 


1.017 


2.837 


3.100 


0.322 


l.E-2 


l.E-2 


7V = 4 


El 


-210.72 


-44.149 


-20.15 


-16.227 


-.463 


0.051 


0.406 


0.468 




E2 


-150.28 


-17.423 


-0.512 


-0.417 


1.357 


1.424 


1.522 


1.571 




Es 


-109.37 


-15.430 


+1.336 


1.355 


+2.979 


3.319 


4.059 


4.087 




<P> 


29.92 


28.76 


1.120 


3.647 


0.397 


0.242 


6.4E-2 


8.7E-3 


N^7 


El 


-210.69 


-44.207 


-20.15 


-16.106 


-.830 


-0.435 


-0.307 


-0.283 




E2 


-150.32 


-17.585 


-0.531 


-0.422 


1.353 


1.407 


1.513 


1.551 




Es 


-109.32 


-15.310 


+ 1.334 


1.353 


+3.019 


3.316 


4.038 


4.055 




<P> 


29.95 


28.92 


1.130 


3.777 


0.721 


8.76E-2 


2.3E-2 


4.7E-3 



Table 3: The same for the levels with {J'^,T) = (2+,0) 





X (MeV) 


10' 


lO'* 


10^ 


10^ 


10^ 


10^ 


N <2 


El 


-17.19 


-16.15 


-14.55 


-12.103 


-11.176 


-11.04 




E2 


0.944 


0.95 


0.952 


0.965 


1.022 


1.079 




<P> 


0.923 


0.95 


1.024 


0.807 


0.136 


2.E-2 


N^3 


El 


-17.213 


-15.665 


-11.142 


-5.953 


-3.950 


-3.80 




E2 


1.029 


1.032 


1.034 


1.047 


1.141 


1.19 




<P> 


1.024 


0.918 


2.842 


1.495 


5.E-2 


8.5E-3 


N 


El 


-17.433 


-15.725 


-8.578 


1.042 


1.082 


1.157 




E2 


1.029 


1.032 


1.034 


1.694 


2.656 


3.038 




<P> 


1.075 


1.095 


5.85 


9.1E-3 


2.9E-2 


2.8E-2 


N^7 


El 


-17.361 


-15.649 


-8.243 


1.042 


1.086 


1.162 




E2 


1.030 


1.032 


1.034 


1.475 


2.524 


2.643 




<P> 


1.062 


1.127 


5.974 


1.2E-2 


3.1E-2 


2.7E-2 



Tabic 4: Experimental results and theoretical estimations of energies of the ground and 
lowest excited states of the ^^C nucleus in MeV 





Exp. 


BFW aa-potential 


MN effective A^A^-potential 


BFW act-potential 




[10] 


[11] 


[6] 


with OPP 




-7.275 


-12.5 


-10.43 


not converged 




0.3796 


-2.7 


0.64 


1.55 




3.0 




5.43 


4.055 




-2.836 


-10.0 


-7.63 


not converged 


n 


3.89 


-3.7 


6.39 


2.61 



